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ABSTRACT 

One of the key aspects of a detector material for space-borne hard X-ray and gamma-ray telescopes is the rate of 
prompt and delayed background events generated inside the material by charged and neutral particles striking 
the detector. These particles are Cosmic Rays, particles trapped in Earth's magnetic field, and secondaries 
from Cosmic Ray interacting with the atmosphere and the spacecraft. Here, we present a preliminary study of 
Cadmium Zinc Telluride (CZT) and its behaviour in space environments. We have used the simulation package 
MGGPOD to estimate the background of the CZT detectors in the proposed Energetic X-ray Imaging Survey 
Telescope (EXIST) for possible orbital parameters. The EXIST mission will make use of ~6 m 2 of >0.5 cm 
thick CZT detectors to record cosmic X-rays in the energy range from 10 keV to 600 keV. The detectors will be 
shielded by a fully or partly active shield. For the specific detector and shielding geometry considered here and 
an orbit with a low (7°) inclination, the background rate is dominated by diffuse extragalactic photons below 
^200 keV. Prompt and delayed hadronic backgrounds grow increasingly important above this energy, becoming 
the main contributors to the total background above ~1 MeV. A fully active shield performs slightly better than 
a half active/half passive shield. 
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1. INTRODUCTION 

With the launch of the Swift satellite in late 2004, CdZnTe (CZT) detectors entered the realm of space based X- 
ray and gamma-ray astronomy. Similar to CZT, CdTe has also been used as a detector material for instruments 
(IBIS and ISGRI) on the INTEGRAL mission prior to the operation of Swift. CZT offers many attributes that 
make it an excellent choice for a hard X-ray survey mission. Its wide band gap allows for operation without 
cryogenic cooling. Its high atomic number (Z = 50) results in high stopping power and a high photo-effect cross 
section. Furthermore, it offers good spatial and energy resolution from 10 keV to 600 keV. 

The sensitivity of hard X-ray telescopes is limited by statistical fluctuations of the background rate, and 
thus the energy dependent background rate is an important property of a detector material. One commonly 
distinguishes between prompt and delayed background events. Prompt background occurs when an incident 
particle immediately creates a signal either by depositing energy or by causing a prompt nuclear reaction. All 
emission occurring within 1 /isec of an interaction is here considered to be prompt. Nuclear reactions can produce 
excited nuclei with long half-lives in the detector and surrounding materials. When these nuclei de-excite, their 
products can deposit energy in the detector producing a delayed signal. This process, referred to as activation, 
has a cumulative effect and the delayed background rates can increase during the lifetime of a mission. 

There are several methods to suppress the background in space-borne telescopes. One is to passively shield 
the detector from incident particles, reducing the flux seen by the detector. While cost-effective and simple to 
design, there are limits to the impact of this technique. Above a certain energy (depending on shield material 
and geometry), particles will still be able to penetrate the shield (shield leakage) and interact with the detector. 
In addition, the shields will become activated from Cosmic Rays and neutrons, producing additional background 
events. In some cases, a better background suppression can be achieved with active shields that generate signals 
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Figure 1. Concept design for the EXIST satellite that can monitor the entire X-ray sky every 95 min orbit. The high 
energy telescopes (HETs) are sensitive from 10 keV - 600 keV. The low energy telescopes (small boxes adjacent to the 
HETs, not labeled) detect photons with energy between 5 and 30 keV. 



when hit. The veto-signals can be used to tag events if one or more particles deposit some of their energy in 
the shield. Furthermore, it is possible to tag "internal background events" generated inside the main detectors 
if high-energy photons leave the detector material and strike the shield. The tagged events can be excluded 
from the data analysis. Active shields require additional electronics and consume power. Some active shields 
add background themselves producing prompt or delayed gamma-ray emission. An optimal shield may use a 
combination of both, active and passive components. 

The Energetic X-ray Imaging Survey Telescope (EXIST) is a proposed satellite-borne hard- X-ray telescope 
that can survey the entire 3-600 keV sky once every 95 min orbit (Grindlay et al. 2006). One design (Figure 1) 
combines a high-energy (HE) CZT detector with a low-energy (LE) Si detector, both being coded mask, wide field 
of view telescopes. The HE telescope will make use of ^6 m 2 of >0.5 cm thick CZT detectors to detect 10-600 
keV X-rays. The LE telescope complements the HE telescope by detecting 3-30 keV X-rays and improving the 
source localization accuracy from 1 arcmin to 10 arcsec. The scientific objective of EXIST is to survey the sky 
for the emission from obscured and non-obscured active galactic nuclei (AGNs), gamma-ray bursts (GRB), and 
galactic black hole systems. A unique feature of the mission with its large field of view would be the capability to 
discover rare transient phenomena, like the catastrophic event of a star plunging into a black hole. EXIST was 
recommended as one of the three high energy missions in the 2000 Decadal Survey Report, and is a candidate 
for the Black Hole Finder Probe (BHFP), one of three Einstein Probe missions recommended in NASAs Beyond 
Einstein Program. 

We are carrying through a program to model the background in balloon and space-borne CdTe and CZT 
detectors. Our primary objective is to optimize the design of EXIST. We plan to verify the modeling by simulating 
the background in the CZT and CdTe detectors of balloon and space borne-missions and by comparing the 
simulated with actually measured background rates and energy spectra. In this contribution, we describe our 
efforts to improve the simulation code and we show first preliminary results from simulating the background of 
the EXIST CZT detectors. 




Figure 2. We have modeled this approximate 
geometry for EXIST, but the actual design 
may be different. The figure shows a side view 
of three high energy sub-telescopes with coded 
mask, CZT detectors, and shield/collimator. 
The shadow pattern produced by the coded 
mask is seen by the CZT detectors. This pat- 
tern can be deconvoluted to reconstruct an X- 
ray image. The active shield/collimator de- 
fines the field-of-view for the CZT detectors 
and reduces the instrumental background . 
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2. SIMULATION DETAILS 

The EXIST HE instrument is made of 19 sub-telescopes (see Figure 1). In rough approximation, we treat each 
sub-telescope as an rectangular box open on one side and closed on 5 sides (see Figure 2). In the following 
we refer to the open side of the collimator/shield assembly as "top" and the closed side as "bottom". Each 
sub-telescope uses a bottom shield and four side slabs for collimation and shielding. The bottom shield is always 
made of CsI(Na). We consider two different designs of the sides, (i) fully active made of CsI(Na), and (ii) half 
active and half passive made of CsI(Na) and Pb-Sn-Cu. In the latter case, the active part surrounds the CZT 
detectors. We use three different veto configurations (no veto, 100 keV threshold, and 1 MeV threshold). The 
bottom shield is active for all simulation runs using veto capabilities. While we assumed here a design based on 
right angles, the actual EXIST design joins shields at angles slightly larger than 90°, resulting in a somewhat 
better photon detection efficiency 

We modeled an approximate geometry for one EXIST HE sub-telescope, as shown in Figure2. Positioned 
directly above the bottom shield is a 1 cm thick region filled with readout electronics, above which lies a 0.5 cm 
thick CZT detector plane, both with an area of 56 x 56 cm 2 . Located ~1.5 m above the CZT is a 0.7 cm-thick 
tungsten coded mask. While the actual mask pixels will have a unit cell size of 2.5 mm x 2.5 mm, we used a 
checkerboard pattern of tiles and gaps, each with an area of 80 x 80 cm 2 , to reduce the simulation time. In 
addition, the coded mask used here provides similar detector response for alternating regions of the field of view 
(FoV). This is in contrast to the actual mask design which is optimized forboth uniform X-ray throughputs in 
the fully coded FoV (22° x 22°) and severe auto-collimation in the partially coded FoV. For the diffuse cosmic 
X-ray background, whose event rates are strongly dependent on FoV, we approximate these effects with a slightly 
different mask model. The FoV up to 22° x 22° contains a checkerboard of 7 mm-thick tungsten mask and gaps 
with a unit cell size of 14 cm x 14 cm. The remainder of the FoV defined by the side shields/collimators (66° 
x 66° FWHM) is covered by continuous 3.5 mm-thick tungsten. This approximates the attenuation effect of the 
partly coded FoV and will be replaced in future simulations by an exact ray-tracing treatment for the radial mask 
holes. We use Al composite bars under the bottom shield to represent support structures for the sub-telescope. 
Finally, we introduce an Al-composite cube to represent the mass of the spacecraft. Table 1 gives the dimensions 
and properties of all components of the mass model. 

We use the MGGPOD software suite (Weidenspointner et al. 2005) to expose the mass model to the dif- 
ferent radiation fields found in a low earth orbit. MGGPOD is comprised of five closely integrated packages: 
MGEANT simulates particle transport and electromagnetic interactions; GCALOR is responsible for handling 
hadron physics; PROMPT handles prompt emission from the de-excitation of nuclei due to spallation from in- 
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[g cm -3 ] 


lxwxt [cm 3 ] 


Detector 


CZT 




5.78 


56 x 56 x 0.5 


CsI(Na) Shields 


Bottom 


CsI(Na) 


4.51 


58 x 58 x 2 


Big Side 


CsI(Na) 


4.51 


58 x 45 x 1 


Small Side 


CsI(Na) 


4.51 


56 x 45 x 1 


CsI(Na) + 


Pb-Sn-Cu Shields 






Big Side Lower 


CsI(Na) 


4.51 


58 x 22.5 x 1 


Big Side Upper 


Pb-Sn-Cu 


10.5 


58 x 22.5 x 0.5 


Small Side Lower 


CsI(Na) 


4.51 


56 x 22.5 x 1 


Small Side Upper 


Pb-Sn-Cu 


10.5 


56 x 22.5 x 0.5 


Other 


Electronics 


Si-compositc 


1.8 


56 x 56 x 1 


Support Bars 


Al-composite 


0.1667 


58 x 1 x 1.5 


Spacecraft 


Al-composite 


0.3675 


124 x 124 x 124 


Mask Tile (1 of 8) 


Tungsten 


19.3 


80 x 80 x 0.7 



Table 1. The densities (g/cm 3 ) and dimensions (length x width x thickness) are given for structures used in the sub- 
telescope mass model 
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65 Cu 














JENDL-3.3 


Cd 


106 Cd 


108 Cd 


110 Cd 


m Cd 


112 Cd 


113 Cd 


114 Cd 


116 Cd 


ENDF/B-VI.8 
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204p b 
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127j 
















JENDL-3.3 


Cs 


133 Cs 
















JENDL-3.3 


Zn 


64 Zn 


66 Zn 


67 Zn 


68 Zn 


64 Zn 


70 Zn 


(Derived) 




JEFF-3.1/A Brond-2.2 



Table 2. Elements and their isotopes for which neutron cross section reformatted and derived (Zn) from the officially 
available data bases. 



cident protons and neutrons, inelastic neutron scattering, and neutron capture; the ORIHET package calculates 
the build up and decay of activity based on nuclide production rates; DECAY then simulates the radioactive 
decays. MGGPOD results match well with the experimental data from the Transient Gamma-Ray Spectrometer 
(TGRS) experiment. 8 

An important part of simulating the delayed background is the consideration of hadronic interactions. The 
GCALOR code consists of 4 packages whose use depends on incident particle and energy. Table 2 summarizes 
these packages and their relevant particles and energies. Of special consideration is the MICAP energy regime, 
as all physics of neutron interactions must be provided via cross-section files. 

Many isotopes of interest had been missing from the original GCALOR/MICAP neutron cross section data 
banks. New cross sections data files were developed by processing data from JENDL-3.3 and ENDF/B-VI.8 data 
sets to make them compatible with the MICAP format. In addition, clement Zn is missing from the JENDL-3.3 



MGEANT 


Particle transport, electromagnetic interactions 


GCALOR 


CALOR package for hadron physics 

1. NMTC (Nuclcon Meson Transport Code) 

nucleons 1 McV to 3.5 GeV, charged pions 1 MeV to 2.5 GeV 

2. Scaling Model 

nucleons and charged pions 3 GeV to 10 GeV 

3. FLUKA 

nucleons and charged pions above 10 GeV and for all energies for particle types 
not implemented in other portions 

4. MICAP 

neutrons 1 • 10" 5 eV to 20 McV 


PROMPT 


Simulates prompt emission from de-excitation of nuclei due to spallation from 
incident protons and neutrons, inelastic scattering and neutron capture 


ORIHET 


Calculates build up and decay of activity based on nuclide production rates 


DECAY 


Simulates radioactive decays 



Table 3. This table describes the five main software components incorporated into MGGPOD. The four software packages 
used by GCALOR are also given. For each GCALOR package, the relevant particles, interactions, and energy regimes 
are given. 

and ENDF/B-VI compilations. Additional neutron cross section files were assembled based on data from the 
JEFF-3.1A, BROND-2.2, and EXFOR data bases. Table 2 lists the isotopes for which cross section files were 
compiled. 

3. BACKGROUND INPUT SPECTRA 

We consider six background components relevant to a 550 km orbit at 7° inclination, other orbital inclinations 
will be considered later. This orbit is realistic for a spacecraft weight of 8500 kg and a launch vehicle of the class 
Delta IV(4050H). Each component is simulated separately, with the final results being combined to form a total 
background. The particles are generated with a random position and trajectory on the surface of a cylinder with 
radius 1500 cm and height 2000 cm, centered about the sub-telescope. Prompt energy deposits in the CZT were 
recorded and filled into histograms, weighted for a 1 sec interval. At the same time, the production of excited 
nuclei was also recorded. Any prompt de-excitations were included in the prompt deposited energy histograms. 
The resulting nuclei production rates were then calculated and evaluated for various orbital histories. The nuclei 
were allowed to de-excite and decay for 30 sec, while the resulting delayed energy depositions in the CZT were 
recorded. 

The simulated background components include diffuse extragalactic X-rays, cosmic ray protons and electrons, 
trapped protons, and atmospheric neutrons. We use the diffuse extragalactic X-ray background as measured by 
the HEAO-1 satellite (Gruber et al. 1999). A flux of atmospheric neutrons (ANs) is produced from cosmic 
ray interactions. We use one of two input spectra, depending on neutron energy: COMPTEL observations 
from Morris (95,98) above 10 MeV, and theoretical calculations from Lingenfelter (1963) below 10 MeV. ESA's 
SPace ENVironment Information System (SPENVIS)* determined the trapped proton peak fluxes incorporated 
in the simulations and CREME96 t computed the time averaged cosmic ray proton (CRP) flux. We assume 
the instruments arc turned off during passages through the South Atlantic Anomaly (SAA), so only delayed 
events are considered for trapped protons. Using the orbital parameters from above, the results from SPENVIS 
show that SAA passages occur typically 3-4 times per day, with an average duration of 10 minutes. For our 
simulations, we assumed SAA passages occur every 6 hours for 10 minutes. Finally, the cosmic electron flux 
from Gchrels (1985) was included with a low-energy cutoff at 2.5 GeV. 

*http:/ /www.spenvis. oma.be 
t https://creme96.nrl.navy.mil/ 
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Figure 3. Total Background 

(prompt+delayed) before and after applying a 
shield- veto whenever an event deposits >100 
keV in any one of the shields (fully active 
configuration). The relative contribution from 
different environmental components are also 
shown (before veto). We assumed an in-orbit 
time of 1 year with activation by primary CR 
protons, albedo neutrons, SAA particles, and 
secondary particles produced by interactions 
of the primary particles with the satellite. 
The figure assumes that 5 min passed since 
the last SAA passage. 




Figure 4. Structural contributions to the de- 
layed background due to one year exposure to 
CR protons and albedo neutrons (after shield- 
veto) . These results are for a fully active shield 
with a 100 keV veto energy threshold. The 
CZT detectors and CsI(Na) shields are sig- 
nificant sources of delayed background events. 
The contribution from the coded mask, elec- 
tronics, and support structures are negligible 
and are not included in this figure. 



4. RESULTS 

The average delayed background rates are expected to increase during the first weeks and months of the mission. 
After roughly half a year, the background rates reach a rather stable mean level with short-time variations owing 
to SAA passages. Figure 3 shows the total (prompt+delayed) background rate in the CZT detectors one year into 
the mission. We assumed a fully active shield and show the total rate before and after applying a veto with an 
energy threshold of 100 keV deposited in any one shield. The vetoes reduce the background rate above 200 keV 
by a factor of 2. The figure also shows the relative contribution from different background components (before 
veto). The dominant sources are diffuse photons and cosmic ray protons. In our model, diffuse photons are most 
important up to ^200 keV. At higher energies, hadronic prompt and delayed events grow increasingly important. 
Above ^2 McV, the CRP prompt events become dominant. The cosmic electron event rate is negligible at all 
energies of interest. The total un-vetoed event rate between 10 and 600 keV is ~0.45 cm' 2 s _1 . 

Figure 4 shows that the delayed background originates mainly in the CZT detectors and in the shield (fully 
active shield, after shield- veto, veto energy threshold of 100 keV). The mask, electronics, and support structures 
contribute little to the delayed signal. The masks are too far away from the CZT detectors. The effects from 
activation of the support structures and electronics is limited due to shielding by the CsI(Na). 

We estimated the delayed signal for different durations in orbit. Figure 5 shows the background rates for 
one hour, one day, six months, one year, and two years of exposure to CRPs, ANs, and their secondaries (fully 
active shields, 100 keV veto energy threshold). One can see the large increase in rate between one month and 
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Figure 5. The delayed background resulting 
from cosmic ray proton and albedo neutron 
interactions are shown for different durations 
in orbit. The background rates increase for 
increasing durations in orbit, but tend to level 
off by six months. Also shown is the delayed 
background due to two years of SAA passages 
5 min after the last SAA passage. 
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Figure 6. This figure shows delayed back- 
ground rates produced by a single SAA pas- 
sage for different "cool down times": just ex- 
iting the SAA (solid thin line), f min (dotted- 
thick line), 5 min (dashed line), 30 min (dot- 
dash line), and 24 hours (dotted-thin line) af- 
ter exiting the SAA. We compare these back- 
ground rates with the long-lived activation 
from the SAA passages accumulating over one 
year in orbit(solid-thick line), shown here 24 
hours after the last passage. These results as- 
sume only passive CsI(Na) shielding. 
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six months, however, the difference between one year and two years is not as dramatic, as nuclei production and 
decay rates tend towards equilibrium. This is expected as there are only a few isotopes with lifetimes longer than 
one year. The same figure shows the delayed event rate from activation by SAA particles and their secondaries 
5 min after the last SAA passage. All three activation components (CRPs, ANs, SAA) are important and 
contribute to the total rate of delayed emission. 

Figure 6 compares the "fast" and "slow" components of the SAA delayed background. Here, fast (FC) refers 
to the events that are products of secondary nuclei with short (<f Day) half-lives, while events from excited 
nuclei with longer half-lives are deemed slow (SC). For one SAA passage, the event rates are shown for various 
cool down times after the passage (CsI(Na) shields with no veto). The cool down times range from immediately 
after exiting the SAA to 24 hours. The long component represents 1 year of SAA exposures. These results show 
that one year into the mission, the long-lived component produces a similar background rate as the short-lived 
component. 

Figure 7 displays the combined prompt and delayed (after one year) background rates using four different 
shielding scenarios; fully active (Csl(Na)) shield with no veto, veto with f 00 keV and f 000 keV energy threshold, 
partly active (Csl(Na)+Pb-Sn-Cu) shield with no veto and veto with 100 keV energy threshold. Delayed CZT 
background rates resulting from CRP, AN, and delayed SAA (last six passages, 5 minutes after last passage) 
events are given. As expected, using a "low" 100 keV shield threshold minimizes the background rate for both 
material configurations. The fully active shield (FAS) provides a lower background rate over the energy range 




Figure 7. The total prompt and de- 
layed background rate is shown for five dif- 
ferent shield/veto-energy-threshold configura- 
tions. FAS stands for "fully active shield" and 
PAS stands for "partly active shield" . The de- 
layed components include contributions from 
CR protons and albedo neutrons after one year 
in orbit, as well as SAA effects from the six 
most recent passages (30 hours) in orbit, 5 min 
after the last SAA passage. 



of interest. For some energies, the detected rate using shield veto is higher than without. This is an effect from 
simulating with different incident photon sets. 

5. CONCLUSIONS 

The preliminary results from our simulations indicate that active shielding is an effective means of reducing the 
background for the EXIST mission. Fully active shields appear to perform better than partly- active/partly- 
passive and fully passive shields. For our detector, mask, and shielding geometry, and assuming a low 7° 
inclination orbit, the EXIST background seems to be dominated by aperture flux (photons entering the detector 
assembly through the optical path), and shield leakage (photons passing through the shield without causing a 
shield- veto) below 200 keV. Above ~1 MeV, the prompt and delayed signals from hadronic background com- 
ponents are the main contributors to the total background rates. Future work will build on these first results, 
incorporating more detailed mass models and environmental conditions. Other detector and shielding configu- 
rations will be evaluated. Another possible design uses more collimation and canted side shields. We plan to 
investigate further the effects from different ratios of active/passive shielding, as well as evaluating the pay-off of 
using thicker shields to reduce the shield leakage. Furthermore, we plan to use the MGGPOD code to simulate 
the background in the CZT/CdTe detectors of the existing missions Swift and INTEGRAL and of past balloon- 
borne experiments. Comparison of simulated with experimentally measured background rates are instrumental 
for checking the simulation tools and the input spectra. 
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